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The  present  analysis  considers  a  three-dimensional  non-isothermal  model  in  a  single  phase  of  a  PEM  fuel 
cell  with  a  flow  field  path  in  the  shape  of  1 , 2, 3, 4, 6,  and  8  concentric  spirals.  The  current  density  contours, 
the  water  content  and  the  entropy  generated  in  all  zones  of  the  fuel  cell  are  predicted.  The  analysis  of  the 
three-dimensional  model  includes  the  gas  flow  channels  in  the  six  geometric  shapes  mentioned  above, 
the  current  collectors,  gas  diffusion  layers,  catalyst  layers  on  both  sides  of  the  model,  anode  and  cathode, 
and  a  proton  exchange  membrane  in  between.  The  energy  equation,  mass  conservation,  and  transport  of 
species  equations  are  solved,  including  source  terms  that  take  into  account  the  electrochemical  effects 
occurring  inside  the  cell.  Also,  the  entropy  generation  equation  is  added  to  the  governing  equations  of  the 
model.  The  results  allow  a  comparison  to  help  to  decide  which  of  the  6  analyzed  configurations  improve 
the  performance  of  the  fuel  cell,  increasing  the  current  density  produced,  reducing  the  pressure  drop 
and  producing  the  most  uniform  current  density.  The  entropy  generation  analysis  reveals  the  effects  that 
cause  the  most  significant  losses  (irreversibilities)  in  the  cell.  The  Bejan  number  and  the  n  number  are 
used  to  compare  the  irreversibilities  produced  by  the  matter  flow  and  by  the  heat  transfer  for  each  one 
of  the  six  models. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Fuel  cells  are  electrochemical  devices  that  convert  the  chemical 
energy  of  a  fuel  directly  to  electrical  energy.  There  is  a  whole  family 
of  fuel  cells  that  can  be  characterized  by  the  electrolyte  used.  All  of 
these  fuel  cells  function  in  the  same  basic  way.  At  the  anode,  a  fuel 
(usually  hydrogen)  is  dissociated,  while  at  the  cathode  oxygen  is 
reduced  to  oxide  species. 

The  proton  exchange  membrane  fuel  cell,  PEMFC,  derives  its 
name  from  the  special  membrane  used  as  the  electrolyte.  The  elec¬ 
trochemical  reactions  that  take  place  inside  the  cell  are: 

H2^  2H+  +  2e“  (1) 

l02  +2H+  +2e~  -*  H20  (2) 

Proton  exchange  membrane  fuel  cells  have  gained  a  lot  of  pop¬ 
ularity  amongst  scientists  due  to  their  operating  conditions,  high 
density  of  energy  generated  and  the  portable  capabilities  of  these 
devices. 
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Fig.  1  shows  the  cross  section  of  a  PEM  fuel  cell.  The  main  com¬ 
ponents  of  a  PEMFC  are  the  current  collectors,  the  gas  diffusion 
layers  (GDL’s),  the  catalysts  layers  and  the  proton  exchange  mem¬ 
brane.  The  main  factors  that  determine  the  fuel  cell  performance 
are  the  operating  conditions,  the  reactant’s  humidity,  the  proper¬ 
ties  of  the  materials  and  the  flow  distribution  of  the  gases  in  the 
MEA  by  means  of  the  geometry  of  the  gas  channels. 

The  geometries  of  the  bipolar  plates  have  been  the  object  of 
many  studies.  The  serpentine  shape  [1,2]  has  proved  to  have  a  good 
performance  as  a  flow  distributor,  with  the  inconvenient  of  produc¬ 
ing  larger  pressure  drops  caused  mainly  by  the  numerous  sudden, 
sharp  turns.  The  parallel  channels  geometry  [3]  is  a  bad  gas  distrib¬ 
utor  causing  larger  concentration  gradients.  The  geometry  called 
interdigitated  [4,5]  is  also  a  deficient  gas  distributor  because  it  pro¬ 
duces  flooding  zones  (although  its  performance  improves  when 
the  operating  pressure  increases.)  These  3  geometries  are  already 
considered  as  conventional  due  to  their  wide  commercial  use. 

Recently,  new  and  non-conventional  flow  configurations  for  the 
flow  field  of  the  PEMFC  have  been  reported  in  the  technical  litera¬ 
ture.  These  results  reveal  the  fact  that  the  fuel  cells  require  a  multi 
objective  optimization  to  make  them  an  option  to  replace  fossil  fuel 
based  technologies.  The  main  objective  behind  those  designs  is  to 
improve  the  fuel  cell  performance,  developing  larger  current  den¬ 
sities  in  smaller  areas.  These  designs,  such  as  the  radial  [6,7],  the 
constructal  theory-based  [8,9],  the  biomimetical  [10],  the  tree  net- 
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Nomenclature 

A  Area  (cm2) 

ak  Water  activity  on  stream  k 

Be  Bejan  number 

Dt  Diffusivity  of  the  i  species  (m2  s-1 ) 

F  Faraday  constant  (96,487  C  mol-1 ) 

fareact  Heat  of  water  formation  ( W  m-3 ) 

I  Current  density  (A  m-2 ) 

j)  Diffusion  mass  flux  of  the  i  species  (kg  m-2  s-1 ) 

jrkef  Reference  current  density  on  stream  k  (A  m-3 ) 
Mi  Molecular  weight  of  the  i  species  (kg  kmol-1 ) 
n  Number  of  turns 

nd  Electrosmotic  drag  coefficient 

P  Total  pressure  (Pa) 

Q  Volumetric  flow  rate  (m3  s-1 ) 

q  Heat  rate  (W) 

R  Universal  gas  constant  (8.314  kj  kmol-1  K-1 ) 

Rk  Exchange  current  density  (A  m-3 ) 

R0hm  Ohmic  resistance  ( £2  m) 

T  Temperature  (I<) 

V  Velocity  (ms-1) 

Voc  Open  circuit  voltage  (V) 

Vcen  Cell  voltage  (V) 

V\/eiec  Electric  power  density  ( W  cm-2 ) 

X/  Molar  fraction  of  the  i  species 

Yi  Mass  fraction  of  the  i  species 

Greek  symbols 

a  Transfer  coefficient 

|3  Permeability  (m2) 

s  Porosity 

c/)k  Potential  field  in  the  k  interface  (V) 

y  Concentration  coefficient 

iq  Overpotential  (V) 

X  Water  content 

p,  Dynamic  viscosity  ( kg  s  m-2 ) 

El  Pi  number 

p  Density  (kg  m-3) 

a  Electric  conductivity  ( £2-1  m-1 ) 

gs  Entropy  generation  (W  K-1 ) 

£  Stoichiometric  flow  rate 

B.  Subscripts/ superscripts 

a  Anode 

act  Activation 

c  Cathode 

cat  Catalyst 

i  H2, 02,H20 

HT  Heat 

H2  Hydrogen 

H20  Water 

k  Anode  or  cathode 

MD  Mechanical  dissipation 

MEA  Membrane  electrode  assembly 

MF  Matter  flow 

mem  Membrane 

02  Oxygen 

ref  Reference 

sat  Saturation 

sol  Solid 

w  Water 

wv  Water  vapor 


Membrane 


Current  Collector 

Anode 

Channel 

GDL  Anode 

Catalyst  Anode 

< — 

<> 

«  4-/%  1  T  ^  /\  /J 

Catalyst  Cathode 


GDL  Cathode 

Cathode 

Channel 

Current  Collector 

z3 

z2 


^0 


Fig.  1.  Cross  section  of  the  fuel  cell.  Regions  where  the  source  terms  are  defined. 

work  [1 1,12],  the  spiral  [13,14],  and  the  fractal  [15]  are  named  after 
the  geometric  shape  or  the  theory  used  to  build  them.  Although 
these  new  designs  present  improvements  they  still  have  deficien¬ 
cies  in  comparison  with  the  conventional  ones.  Amongst  these 
deficiencies  are  the  poor  active  area  utilization  which  leads  to  low 
current  density  production;  and  high  ohmic  losses  due  to  a  defi¬ 
cient  contact  area  between  the  current  collectors  and  the  electrode. 
These  deficiencies  are  inherent  to  the  structural  design. 

All  real  processes  present  irreversibilities  that  can  be  associated 
to  friction,  mass  transfer,  heat  transfer,  electrochemical  reactions, 
etc.,  resulting  in  a  non  desirable  loss  or  degradation  of  energy  qual¬ 
ity.  This  loss  of  efficiency  is  related  to  the  entropy  generation.  The 
generation  of  entropy  is  a  synonym  for  the  potential  destruction 
of  system  energy;  its  minimization  has  been  used  as  the  optimal 
design  criteria  for  thermal  systems.  The  mechanisms  of  entropy 
generation  must  be  understood  [16,17]  by  analyzing  the  process 
itself  and  its  thermal,  hydrodynamic  and  geometric  characteristics 
in  order  to  obtain  an  optimal  configuration  with  a  minimum  loss 
of  available  energy.  Damian  et  al.  [18]  presented  the  entropy  gen¬ 
eration  analysis  of  a  PEMFC  with  a  biomimetical  flow  field.  Their 
analysis  consists  of  not  only  the  calculation  of  the  total  entropy  gen¬ 
erated,  but  also  an  analysis  of  the  impact  of  the  entropy  generation 
on  each  of  the  different  processes  involved  on  the  cell  performance. 
Their  results  indicate  that  the  flow  and  thermal  effects  are  smaller 
than  the  entropy  generated  by  the  transport  of  species. 

The  present  work  shows  the  results  for  the  overall  performance 
and  entropy  generation  of  a  PEMFC  with  multiple  spirals  as  flow 
fields.  The  results  of  the  entropy  analysis  are  compared  with  the 
performance  analysis  to  have  a  complete  description  of  the  cell 
performance  in  the  hope  of  finding  the  most  efficient  geometric 
shape. 

2.  Model  development 

The  present  work  is  a  numerical  analysis  of  six  fuel  cells  with 
n  =  1,  2,  3,  4,  6  and  8  flow  channels  of  constant  width,  see  Fig.  2. 
The  geometries  are  generated  using  concentric  Archimedes’  spi¬ 
rals,  r=a0,  with  a  2 7 r/n  radian  phase  angle  between  each  couple  of 
spirals.  This  allows  generating  two  channels  of  constant  width. 

A  computational  three-dimensional  model  is  built  for  each  of 
the  6  geometries  considered.  The  models  are  used  to  analyze  the 
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1  Channel 


2  Channels 


3  Channels 


4  Channels 


6  Channels  8  Channels 

Fig.  2.  Geometry  of  the  fuel  cells  analyzed. 


electrochemical  reactions  and  the  transport  phenomena  occurring 
at  the  catalyst  layers  (at  both  the  anode  and  the  cathode  side).  In 
other  areas,  such  as  the  gas  diffusion  layer,  no  chemical  reactions 
occur. 

The  geometric  dimensions  are  listed  in  Table  1.  These  dimen¬ 
sions  are  similar  to  the  dimensions  used  by  Shimpalee  and  Dutta 
[  1  ].  The  current  collectors,  GDL,  catalysts  layers,  and  the  membrane 
have  the  same  width  as  the  cell. 

The  gas  flow  channels  and  the  current  collectors  are  the  zones 
of  the  cell  that  are  adapted  to  the  new  geometry,  see  Fig.  2.  The  rest 


of  the  fuel  cell  zones:  GDL’s,  catalyst  layers,  anode  and  cathode,  and 
also  the  membrane,  are  not  affected  by  this  geometry  change.  The 
number  of  turns  varies  in  such  a  way  that  the  reaction  area  is  the 
same  for  all  the  models. 

Fig.  3  shows  the  three-dimensional  model  developed  for  the  8 
channel  fuel  cell,  where  the  gas  reactant  inlets  are  located  at  the 
center  of  the  spirals.  The  outlets  are  located  on  the  outer  part  of  the 
spirals.  It  is  important  to  note  that  the  gas  entrance  is  not  at  the 
same  level  as  the  cell.  A  straight  circular  conduct  is  added  to  con¬ 
sider  the  acrylic  base  upon  which  the  fuel  cell  is  mounted.  The  fuel 
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Table  1  4.  Governing  equations 

Model  dimensions. 


Description 

Dimension 

Current  collector  height 

1.50  x  10-3  m 

Channels  height 

1.00  x  10-3  m 

Channels  width 

0.80  x  10~3  m 

GDL  height 

0.50  x  10-3  m 

Catalyst  width 

0.03  x  10~3  m 

Cell  width 

1.60  x  10~3  m 

Inlet  radius 

0.53  x  10-3  m 

Inlet  channel  length 

3.50  x  10-3  m 

Number  of  turns  (1 -channel) 

13.75 

Number  of  turns  (2-channels) 

6.375 

Number  of  turns  (3-channels) 

4.666 

Number  of  turns  (4-channels) 

3.187 

Number  of  turns  (6-channels) 

2.333 

Number  of  turns  (8-channels) 

1.718 

Reaction  area 

16.6  x  10“4  m2 

Fig.  3.  Boundary  conditions  for  the  3-dimensional  model  of  the  8-channel  fuel  cell. 


cell  cross  section  area  defined  previously  in  Fig.  1  can  be  appreciated 
as  the  spiral  branches  out. 


3.  Model  assumptions 

Steady  state,  laminar  flow,  non-isothermal,  single  phase, 
isotropic  materials,  and  mass  transport  along  three  dimensions 
are  assumed  in  the  analysis.  The  open  circuit  voltage  at  the  cell 
temperature  is  V0c  =  1.0904  V. 


The  governing  equations  including  the  source  terms  are  listed  in 
Table  2.  The  added  source  terms  due  to  the  electrochemical  effects 
are  listed  in  Table  3.  The  ‘z’  indicators  in  Fig.  1  show  the  zones  where 
the  source  terms  are  calculated. 

The  mass  diffusion  flux,  Jy,  in  the  species  transport  equation  is 
defined  as 


r  r»  dml 

Ji'l  =  -pD,9^ 


where  Dj  [19]  is  defined  as, 
■  Y  \  1-5 


D  i  = 


(17) 


(18) 


Here  p0  =  1 01 325  N  rrr2  and  T0  =  300  K.  D*  already  considers  the 
porosity  effects.  The  electrosmotic  drag  coefficient  is  calculated 
using  the  following  equation, 


ad  = 


2.5A 

22“ 


(19) 


while  the  water  diffusion  coefficient  at  the  membrane,  Dw,  is  that 
proposed  by  Nguyen  and  White  [20],  i.e. 


Dw  =  5.5  x  10  1 1  exp 


2416 


(20) 


Table  3 

Source  terms. 


Source  terms  and  location  of  application  (see  Fig.  1) 
Sm  =  Sh2  +  Saw  atZ  —  Z3 
Sm  =  So  +  Sc  w  atZ  =  z2 


(10) 


=  _/xu  x 

;  *  1  t 

bpy  =~~rT  }  3t2i<Z<Z4 


Snz  — 


H  w 


Sh  =  I2Rohm  +  h  react  +  ^act^a; 


at  all  zones 


Sh2  = 

Saw  = 


Mh2 


H2°  , 


F 

Mo7 

So2=-^Rc 

_  Mh2o 

JCW  —  2  p 


at  z  =  z3 
at  z  =  z3 
at  z  =  z2 
at  z  =  z2 


(11) 

(12) 

(13) 

(14) 

(15) 

(16) 


Table  2 

Governing  equations. 
Governing  equations 


Continuity 


Momentum  transport 


Energy 

Hydrogen  transport  (anode  side) 
Water  transport 
(Anode  side) 

Oxygen  transport  cathode  side) 
Water  transport  (cathode  side) 


Mathematical  expressions 
d(pu)  d(pv)  d(pw)  _ 

rv  "T~  rv  '  rv  *^IT1 

ox  dy  oz 

d(pu)  d(pu)  d(pu )  dP  d  (  du\  d  {  du\  d  (  du\ 

u^r+v^+w^r  =  -Yx  +  irA^)  +  Ty^iry)  +  ¥z(^Yz) 

d{pv)  d{pv)  a (pv)  ap  a  /  dv\  a  (  av\  a  (  av\ 

uY^+v^  +  wYar=-Y  +  Yx(,1Yx)  +  Yy('xYy)  +  Yz{llYz) 

d(pw )  3 (pw)  3 (pw)  3 P  3  /  3w\  3/3 w\  3  /  3w 
u^r +VYF~ +  w^r~  =  -&  +  &  ("to ) +  Yy  (»-&)  +  Yz  (* & 
a  (per)  a  (per)  d(pcr)  a  f, ar \  a  (,  ar\  a  (.a t\  _ 

u~lx~  +v~ly~  +W~lz~  =  3x  (fc3x  )  +  dy  \  3y )  +  3z  \  Y )  +Sh 

3  (pmH2)  3[pmH2)  3(yX)H2)  9(/yjH2)  9(jz>h2) 

u  — - L  +v— i - L  +  w— - -  =  — S - L  +  -A - L  +  -A - L  +s 

dx  dy  3 z  3x  3y  3z 


dx 
d(pm 


-  +v 


d(pma 


-  +  w 


3z 

d(pmaw)  9  (Jx,aw) 


d(Jy, 


dy 

d(Jz, 


dx  dy  dz  dx 

dfpmo2)  d  (pmo2)  3  (pm02)  3  (jx,o2 

— - -  +  y — - —  +  w — - —  =  — - 

3x  dy  dz  dx 

)  t  A  d(pmcw)  3( _  _ 

dx  +  3y 


dx 
3 (pma 
dx 


-  +v 


y  dz  dx 

d(pmcw)  d(pmcw )  3(JX,CW)  3C/y,cw) 

3y  3z 


3y  3z 

)  3(iy.o2) 

+  — ^ ~  + 


I 

■  +  Saw 
d  (jz,  o2  ) 


dy 

d  (Jz,  cw) 


3z 


3z 


+  Scv 


(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 
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The  saturation  pressure  is  calculated  with  the  correlation  pro¬ 
posed  by  Springer  [21  ],  in  atm,  as, 


The  total  entropy  production  is  obtained  through  integration  of 
the  local  entropy  generation  (cr)  in  each  computational  domain: 


logioPsat  =  -2.1794  +  0.02953  (T  -  273.17)  -  9.1837 
X  10“5(T  -  273. 17)2  +  1.4454  x  10“7(T  -  273. 17)3  (21) 


The  membrane  phase  electric  conductivity,  crmem,  is  calculated 
using  the  equation  also  proposed  by  Springer  [21  ],  namely, 


Omem  —  £(0.514A 


0.326)  exp 


,26s  (a® 


(22) 


while  the  water  content  is  calculated  with  the  correlation  used  by 
Springer  [21], 


f  0.043  +  17.18cza  -  39.85q2  +  36.0a;*  aa  <  1 
|  14.0  +  1.4(aa  -  1)  aa  >  1 


(23) 


here  aa  is  the  water  activity  defined  as 


Qa  = 


Pwv 

Psat 


(24) 


Two  potential  equations  are  solved  in  the  model;  one  poten¬ 
tial  equation  accounts  for  the  electron  transport  through  the  solid 
materials  and  the  other  represents  the  protonic  transport  of  H+.  The 
potential  equations  are  defined  as 


^  (tfsolV0sol)  +  Psol  —  0 


(25) 


crsTi=  f  crfdV  (34) 

Jv 

An  alternative  parameter  to  represent  the  distribution  of  irre¬ 
versibilities  is  the  Bejan  number  (Be)  defined  as  [17]: 

Be  =  — ^  (35) 

crs 

When  the  entropy  generation  is  due  only  to  mechanical  dis¬ 
sipation  and  heat  transfer,  this  dimensionless  number  indicates 
that  heat  transfer  irreversibilities  are  dominant  when  Be  » 1/2. 
However,  in  fuel  cells  the  entropy  generation  depends  mainly  on 
the  matter  diffusion.  For  this  reason,  the  parameter  proposed  by 
Damian  et  al.  [18]  is  used: 


crs 


(36) 


This  parameter  is  the  ratio  of  entropy  generation  due  to  the 
matter  flow  to  the  total  entropy.  This  number  implies  that  irre¬ 
versibilities  due  to  matter  flow  are  insignificant  when  n  « 1/3. 


5.  Boundary  conditions 

The  boundary  conditions  assumed  for  the  computational 
domain  are  shown  in  Fig.  4.  At  the  inlet  of  the  gas  flow  channels, 
the  boundary  condition  is  given  by  the  stoichiometric  flow  rate: 


V  (CTmemV0  mem)  +  Pmem  —  0 


(26) 


The  transfer  current  or  source  terms  are  nonzero  only  inside 
the  catalyst  layers  and  are  determined  as  Psol  =  -Ra  for  the  solid 
phase  at  the  anode  side  and  Rso\  =  +RC  at  the  cathode  side.  For  the 
membrane  phase,  Rmem  =  +Pa  at  the  anode  side  and  Fmem  =  -Pc  at 
the  cathode  side. 

The  source  terms  in  Eqs.  (26)  and  (27)  are  calculated  via  the 
Butler-Vollmer  equation,  namely, 


(  [H2] 

( <XaFr]a  \ 

/acFpaV 

V  [W2]ref  / 

exp  | 

V  RT  J 

1  6XP  \  -RT  )\ 

f  [02]  Vc 

/’acFrjc\ 

(  a3Fr]c  \ ' 

( [  02  ]  ref  / 

exp  ( 

v  -RT  ) 

”exp(  RT  JJ 

(27) 

(28) 


The  overpotentials  at  the  anode  and  cathode  sides  are  found 
from 


IJa  —  0sol  —  0mem  (29) 

*lc  —  0  sol  —  0mem  —  Voc  (30) 

while  the  open  circuit  reference  voltage  [22]  is  given  by 

Voc  =  0.0025T  +  0.2329  (31 ) 

where  T  represents  the  operation  cell  temperature. 

The  total  entropy  generation  is  calculated  as  follows  [23]: 

jv  o 

as=qv(l)-^kv(^)-lpvVu-Pu:Vu  (32) 

k= l 


From  Eq.  (32)  it  is  clear  that  three  different  effects  contribute 
to  the  rate  of  entropy  production,  namely:  the  entropy  production 
related  to  heat  transfer,  the  second  term  to  the  matter  flow  and 
the  third  and  fourth  to  mechanical  dissipation.  Thus,  the  entropy 
production  can  be  expressed  as: 

°’S  =  °’HT  +  °’MF  +  crMD  (33) 


£a  —  ^H2,inQa,in 


Pa  2F 
RT  /refAMEA 


£c  —  ^O2,in0c,in 


Pc  4F 
RT  /refAMEA 


(38) 

(39) 


A  zero-gauge  pressure  condition  is  defined  at  the  outlets  to 
simulate  a  discharge  to  the  atmosphere.  The  external  faces  of  the 
current  collectors  are  defined  as  a  wall  with  the  condition  of  a 
potentiostatic  cell,  i.e.  the  electric  potential  at  the  anode  is  0sol  =  0 
and  at  the  cathode  0soi  =  Vce\\. 
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Table  4 

Operating  conditions. 


Parameter 

Value 

Pressure 

202650  Pa 

Cell  temperature 

343  K 

Humidity  at  the  anode  side 

100% 

Humidity  at  the  cathode  side 

100% 

Anode  stoichiometric  flow  rate 

1.5 

Cathode  stoichiometric  flow  rate 

2 

Table  5 

Electrochemical  parameters. 

Property 

Value 

Current  collector  conductivity  [25] 

orsoi  =  1.25  x  105  J2-1  m-1 

GDL  and  catalyst  conductivity  [25] 

O'  cat  =  O'  GDL  =  53  ITT1 

GDL  porosity  [24] 

£gdl  =  0.4 

Catalyst  porosity  [24] 

Scat  =  0.1 12 

Concentration  exp.  (anode)  [24] 

Y  a  =  0.5 

Concentration  exp.  (cathode)  [24] 

Kc  =  l 

Transfer  coefficient  (anode)  [24] 

Q'a  =  2 

Transfer  coefficient  (cathode)  [24] 

ac  =  2 

Ref.  current  density  (anode)  [24] 

jf  =  1  x  109A  m“3 

Ref.  current  density  (cathode)  [24] 

jf  =  3  x  105A  m-3 

Hydrogen  reference  diffusivity  [24] 

D°  =  1.1  x  10“4m2s-1 

H2 

Oxygen  reference  diffusivity  [24] 

D°  =  3.2  x  10~5  m2  s-1 
°2 

Water  reference  diffusivity  [24] 

D°  o  =  7.35  x  10-5  m2  s-1 

Nitrogen  reference  diffusivity  [12] 

D°2  =  8  x  10-5  m2  s_1 

Permeability  [12] 

/3  =  2  x  10-10  m2 

The  terminal  walls  are  defined  as  adiabatic,  while  the  side  walls 
are  defined  at  the  cell  temperature. 

6.  Operating  conditions 

The  model  was  analyzed  under  the  operating  conditions  listed 
in  Table  4. 

The  electrochemical  parameters  employed  in  the  numerical 
model  are  listed  in  Table  5. 

7.  Numerical  procedure 

The  mesh  quality  was  verified  and  refined  in  order  to  avoid  hav¬ 
ing  a  model  with  highly  skewed  elements.  Fig.  5  shows  the  mesh 


quality  distribution  for  the  8-channel  model,  the  one  with  the  larger 
number  of  highly  skewed  elements. 

A  mesh-independence  analysis  was  conducted  varying  the  total 
number  of  elements.  The  grid-independence  was  verified  via  the 
continuity  equation’s  solution.  The  final  mesh  for  the  6  models  had 
an  average  of  1  million  elements. 

A  finite  volume  code  is  used  to  solve  the  coupled  equations  of 
mass  conservation,  mass  diffusion  and  momentum  conservation. 
The  code  includes  a  module  that  takes  into  account  the  electro¬ 
chemical  effects  of  the  source  terms  associated  with  each  zone.  The 
numerical  procedure  is  as  follows:  the  Navier-Stokes  equations  are 
solved  first,  next  a  correction  of  pressure  is  performed  in  order  to 
do  a  correction  in  the  mass  balance,  and  then  the  mass  transfer 
equations  are  solved  with  the  new  data  obtained  from  the  correc¬ 
tion  of  the  mass  balance.  Once  the  properties  have  been  updated 
the  energy  equation  is  solved. 

8.  Results 

8.1.  Validation  of  the  numerical  model 

A  validation  of  the  numerical  model  was  conducted  by  repro¬ 
ducing  the  experimental  performance  curve  reported  in  [24].  Fig.  6 


Quality 


Fig.  5.  Cell  quality  histogram.  Equilateral  volume  deviation  method. 
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shows  the  comparison  between  the  experimental  and  the  numeri¬ 
cal  simulation  using  the  same  operating  conditions  and  dimensions 
presented  by  Um  and  Wang  [24].  The  comparison  shows  a  very 
good  agreement  between  numerical  and  experimental  data.  The 
present  validation  ensures  that  the  model  described  so  far  works 
as  a  PEM  fuel  cell  performance  simulator. 

8.2.  Polarization  and  power  curves 

The  polarization  and  power  curves  for  the  6  fuel  cell  models 
analyzed  are  presented  in  Figs.  7  and  8,  respectively.  The  polariza¬ 
tion  curves  allow  comparison  of  the  overall  performance  of  the  cell 
under  different  operating  conditions. 

The  activation  losses,  which  occur  at  very  low  current  densi¬ 
ties  (0-0.1  A  cm-2),  were  very  similar  for  all  six  models,  see  Fig.  7, 
because  all  the  cell  models  need  almost  the  same  amount  of  energy 
to  begin  the  reaction. 

The  main  difference  in  terms  of  fuel  cell  performance,  amongst 
the  6  geometries  proposed,  is  the  behavior  of  the  ohmic  losses. 
The  ohmic  losses  effect  can  be  easily  observed  by  comparing  the 
polarization  curves  slope.  A  larger  slope  value  means  a  larger  ohmic 
losses  effect.  This  kind  of  losses  appears  at  current  densities  larger 
than  0.1  A  cm-2.  The  ohmic  losses  effect  is  most  significant  for  the 
8-channel  model  having  a  poor  performance.  The  4-channel  and 
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Fig.  7.  Polarization  curves. 
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Fig.  9.  Pressure  drop  for  the  cathode  side. 


the  3-channel  models  have  the  best  performance.  The  2-channel 
and  6-channel  models  have  a  similar  behavior,  although  not  the 
best. 

Finally,  a  comparison  of  the  power  curves,  see  Fig.  8,  supports 
the  previous  discussion.  The  4-channel  fuel  cell  is  the  one  with  the 
largest  power  developed,  while  the  8-channel  fuel  cell  has  the  worst 
performance,  having  a  maximum  power  difference  of  15%  com¬ 
pared  to  the  maximum  power  developed  for  the  4-channel  fuel 
cell. 

8.3.  Pressure  drop 

Fig.  9  shows  the  pressure  drop  on  the  cathode  side  for  each  of 
the  6  models.  As  shown  in  Fig.  9,  the  pressure  drop  increases  as  the 
number  of  spirals  decreases.  Thus,  the  1 -spiral  fuel  cell  produces 
the  largest  pressure  drop  because  the  gases  have  to  go  through  a 
larger  channel  compared  with  the  channel  length  for  the  rest  of  the 
models.  The  cell  voltage  has  a  minimum  effect  on  the  pressure  drop 
having  a  maximum  difference  of  10%  between  the  maximum  and 
minimum  values. 

8.4.  Current  density  contours 

For  a  fuel  cell  it  is  important  to  generate  large  values  of  electrical 
current,  but  how  this  is  done  is  also  important.  Fig.  10  compares 
the  current  density  distribution  for  the  6  fuel  cell  models.  It  is  clear 
that  the  4-channel  fuel  cell  generates  the  electrical  current  in  a 
uniform  way  throughout  of  the  cell,  while  the  8-channel  does  it  in 
an  irregular  way.  The  8-channel  fuel  cell  model  produces  a  larger 
electrical  current  near  the  spiral  center.  However,  at  the  outer  zones 
the  current  produced  is  almost  null  because  the  reactant  gases  have 
been  depleted  when  they  reach  this  zone.  The  rest  of  the  models 
have  a  similar  distribution  to  that  of  the  4-channel  model  except 
for  a  lower  current  density  value. 

8.5.  Distribution  of  gases 

Although  current  density  is  a  good  parameter  to  compare  cell 
performance,  it  is  not  the  only  parameter  that  should  be  taken  into 
account.  The  distribution  of  gases  through  the  cell  is  another  factor 
used  to  examine  the  performance  because  it  affects  the  current 
density  produced,  see  Fig.  1 1 .  A  uniform  gas  distribution  will  extend 
the  lifetime  of  the  membrane  and  will  enhance  the  performance  of 
the  cell.  The  water  produced  by  the  3-channel  fuel  cell  accumulates 
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Fig.  10.  Current  density  produced  (Am-2)  by  the  six  geometries  at  Vcen  =  0.6  V. 


close  to  the  outlets,  which  makes  it  easy  to  remove  it.  However,  it 
also  shows  an  irregular  water  production.  The  8-channel  fuel  cell 
shows  a  more  regular  water  production,  which  will  help  to  extend 
the  PEM  lifetime. 

8.6.  Entropy  generation  by  heat 

Fig.  12  shows  the  entropy  generated  by  each  of  the  fuel  cells 
geometries.  Clearly,  the  entropy  generated  by  heat  increases  as  the 
current  density  increases.  The  heat  is  produced  either  by  the  Joule 
effect  or  the  chemical  reactions  that  take  place  in  the  cell.  The  model 
that  generates  the  least  amount  of  entropy,  for  a  specific  current 
density,  is  the  4-channel  fuel  cell.  On  the  other  hand,  the  fuel  cell 


with  1  spiral  generates  the  largest  amount  of  entropy  by  heat.  The 
rest  of  the  geometries  produce  almost  the  same  amount  of  entropy 
by  heat. 

The  4-channel  fuel  cell  also  produced  the  maximum  entropy  by 
heat  at  higher  current  densities.  Under  this  condition,  the  number 
of  reactions  taking  place  inside  the  cell  is  large;  this  increases  the 
heat  within  the  cell  and  also  the  entropy  produced  by  heat. 

8. 7.  Entropy  generation  by  flow  (viscous  effects ) 

The  entropy  generated  due  to  the  viscous  effects  is  shown  in 
Fig.  1 3.  The  amount  of  entropy  produced  has  no  significant  variation 
compared  to  the  current  density.  A  comparison  between  heat  and 


[W/K]  a  *  [W(K] 


D.  Juarez-Robles  et  al.  /  Journal  of  Power  Sources  196  (201 1 )  8019-8030 


8027 


8.00e-01 

7.60e-01 

7.20e-01 

6.80e-01 

6.40e-01 

6.00e-01 

5.60e-01 

5.20e-01 

4.80e-01 

4.40e-01 

4.00e-01 

3.60e-01 

3.20e-01 

2.80e-01 

2.40e-01 

2.00e-01 

1.60e-01 

1.20e-01 

8.00e-02 

4.00e-02 

0.00e+00 


Fig.  11.  Mass  fraction  of  water  at  the  cathode  catalyst. 


Fig.  12.  Entropy  generation  by  heat. 
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Fig.  13.  Entropy  generation  by  flow. 


flow  entropy  produced  shows  that  the  viscous  effects  are  almost 
negligible. 

The  3-channel  and  4-channel  fuel  cell  produced  the  largest 
quantity  of  entropy  due  to  the  viscous  effects.  However,  this  param¬ 
eter  can  be  considered  negligible  for  two  reasons;  first,  because 
of  its  small  magnitude,  and  second,  because  the  two  geometries 
are  the  ones  that  produce  the  best  cell  performance  (based  on  the 
polarization  and  the  power  curves.) 

8.8.  Entropy  generation  by  matter  flow 

Figs.  14-17  show  the  entropy  generated  by  the  diffusion  through 
the  cell  of  the  hydrogen,  water,  nitrogen  and  oxygen,  respectively. 
The  entropy  generation  by  matter  flow  depends  mainly  on  the  gases 
concentration  and  its  capability  to  diffuse  into  the  porous  zone. 
Thus,  if  the  quantity  of  entropy  generated  is  large,  it  means  that 
there  is  a  large  variation  of  the  gas  concentration  from  one  zone  to 
other. 

The  entropy  generated  by  the  hydrogen  diffusion,  see  Fig.  14, 
shows  that  at  lower  current  density  (voltages  close  to  the  open  cir¬ 
cuit  voltage)  the  entropy  generated  is  larger.  This  is  due  to  the  fact 


MF,02  °  MF,N2  ^  °  MF.H20  ^ 
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Fig.  15.  Entropy  generation  by  the  H20  diffusion. 
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Fig.  16.  Entropy  generation  by  the  N2  diffusion. 


Fig.  17.  Entropy  generation  by  the  02  diffusion. 


that  at  lower  current  densities  the  reactions  occurring  in  the  cell 
are  spontaneous  and  non-uniform.  Under  this  condition,  the  hydro¬ 
gen  enters  the  cell  but  does  not  react,  and  its  chemical  potential  is 
wasted. 

However,  as  the  current  density  increases,  the  entropy  gener¬ 
ated  reaches  a  minimum  point  to  increase  again  its  value  as  the  cell 
voltage  is  close  to  zero.  The  six  geometries  again  produce  almost  the 
same  amount  of  entropy,  within  the  range  of  0-1 A  cm-2,  with  the 
exception  of  the  8-channel  fuel  cell  which  shows  an  erratic  behav¬ 
ior  (due  to  the  non-uniform  current  density  production.)  The  main 
difference  appears  at  larger  current  densities  where  the  4-channel 
fuel  cell  produces  the  lowest  amount  of  entropy  while  the  1  channel 
fuel  cell  produces  the  largest. 

The  entropy  produced  by  the  water  diffusion,  see  Fig.  15,  has  an 
inverse  behavior  than  that  of  the  hydrogen  diffusion  with  a  magni¬ 
tude  impact  10  times  larger.  The  difference  between  the  water  and 
the  hydrogen  entropy  production  can  be  found  by  looking  at  how 
the  cell  works.  At  larger  cell  voltages,  the  water  production  is  small 
because  of  the  small  current  density  required.  As  the  cell  voltage 
decreases,  the  current  density  produced  increases  as  does  the  water 
produced.  This  amount  of  water  produced  should  be  removed  from 
the  cell  through  the  channels.  If  the  geometry  of  the  channels  does 
not  help  this  draining  effect,  the  cell  performance  diminishes.  For 
the  8-channel  and  the  6-channel  fuel  cells,  the  entropy  generated 
due  to  the  water  diffusion  is  large  because  the  channels  cannot 
remove  the  water  surplus.  On  other  hand,  the  1,  2  and  3-channel 
fuel  cells  have  no  problem  with  this  effect  and  they  produce  the 
lowest  entropy  production  by  water  diffusion. 

Nitrogen  is  a  particular  case  on  the  fuel  cells  because  it  enters 
the  cell  joined  to  the  oxygen  as  air,  but  it  does  not  react.  It  can  be 
considered  as  a  cooling  fluid  because  it  takes  heat  from  the  cell 
before  leaving  it.  However,  its  entropy  contribution,  see  Fig.  16,  is 
significant  although  the  mass  of  nitrogen  that  enters  the  cell  is  the 
same  as  that  observed  when  it  exits.  The  nitrogen  concentration 
variation  is  caused  by  the  variation  of  the  rest  of  the  species.  The 
geometries  that  produce  the  larger  amount  of  entropy  by  nitrogen 
diffusion  are  the  6  and  8-channel  fuel  cells.  The  1  and  2-channel 
fuel  cells  produce  the  lowest  amount  of  entropy  due  to  nitrogen 
flow. 

From  all  the  species  that  affect  cell  performance,  oxygen  diffu¬ 
sion  is  the  one  that  produces  the  largest  entropy,  see  Fig.  17.  Once 
again  the  6  and  8-channel  fuel  cells  have  the  larger  entropy  pro¬ 
duction  than  the  rest  of  the  geometries.  Contrary  to  the  effect  of 
the  hydrogen  diffusion,  at  large  cell  voltages  the  entropy  produc¬ 
tion  by  oxygen  diffusion  is  negligible  because  it  does  not  have  to 
react  to  accomplish  the  current  density  demanded.  However,  at  low 
cell  voltage,  the  current  density  produced  increases  considerably, 
producing  large  oxygen  concentrations  variations,  and  therefore 
increasing  the  entropy  produced. 

The  entropy  production  due  to  matter  flow  is  the  sum  of  the 
entropy  contribution  of  H2,  H20,  N2  and  02.  The  total  entropy  pro¬ 
duced  is  the  sum  of  the  contributions  by  means  of  heat,  flow,  and 
matter  flow.  The  results  are  shown  in  Fig.  18  and  they  lead  to  the 
same  conclusion  about  the  best  geometries  (1  and  2  spirals)  and 
the  worst  geometries  (6  and  8  spirals)  with  respect  to  the  entropy 
generation.  Although,  the  performance  of  the  3  and  4-channel  fuel 
cell  is  not  bad  at  all,  developing  almost  the  same  amount  of  entropy 
than  that  of  the  1  and  2  channel  geometries. 

So  far,  the  effects  and  causes  of  each  of  the  processes  (heat,  flow 
and  matter  flow)  that  contribute  to  the  net  entropy  generation  have 
been  compared.  However,  it  is  important  to  decide  which  one  of 
them  is  the  most  significant. 

The  Bejan  number  is  used  to  compare  the  effect  of  the  entropy 
production  by  heat  on  the  net  entropy,  see  Fig.  19.  Bejan  numbers 
lower  than  0.5  mean  that  the  heat  effect  is  not  significant  on  the 
phenomena.  The  maximum  Bejan  number  occurs  at  higher  current 
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Fig.  18.  Total  entropy  production. 


densities  and  it  decreases  as  long  as  the  current  density  decreases. 
In  this  analysis,  the  maximum  Bejan  number  is  0.33  for  the  1  -spiral 
fuel  cell  at  a  cell  voltage  of  0.2  Volts  (JAve  =  1.42  A  cm-2).  Although, 
the  maximum  value  obtained  for  the  Bejan  number  represents  33% 
of  the  net  entropy  produced,  it  does  not  represent  the  main  effect  on 
the  cell  performance.  So,  instead  of  analyzing  the  cell  performance 
using  the  Bejan  number,  it  is  better  to  analyze  it  with  the  entropy 
ratio,  n,  proposed  by  Damian  et  al.  [18],  see  Fig.  20. 

The  n  number  represents  the  ratio  of  the  entropy  generated 
by  matter  flow  on  the  total  entropy  generated.  Lower  n  numbers 
represent  a  more  efficient  fuel  cell.  Fig.  19  shows  that  the  1  and 
2-channel  fuel  cells  are  the  best  geometries  in  terms  of  minimizing 
the  entropy  production  (while  the  6  and  8-channels  are  the  worst.) 

An  analysis  using  the  Bejan  number  will  lead  to  conclude  that 
the  6  and  8  channels  are  the  best  geometries  because  they  produced 
the  lowest  entropy  generation.  Flowever,  the  analysis  using  the  n 
number  will  lead  to  conclude  that  the  1  and  2-channel  fuel  cells 
are  the  best  geometries.  The  way  to  decide  which  one  of  them  is 
more  appropriate  depends  on  deciding  what  is  the  most  important 
effect  in  the  fuel  cell,  heat  production  or  gas  diffusion.  The  main 
function  of  a  fuel  cell  is  to  produce  current  density.  Current  density 


Fig.  19.  Bejan  number. 


Fig.  20.  Pi  number. 


is  related  to  the  number  and  distribution  of  the  reactions  occurring 
in  the  cell.  The  distribution  of  reactions  depends  mainly  on  the  gases 
distribution;  thus,  the  n  number  is  more  suitable  to  characterize 
the  cell  performance. 


9.  Conclusions 

The  present  work  compares  the  performance  of  a  fuel  cell  using 
six  different  geometries,  with  the  Archimedes  spiral  as  a  base.  The 
performance  analysis  is  based  on  the  current  density  produced  and 
its  distribution  along  the  cell.  Also,  an  entropy  analysis  was  con¬ 
ducted  in  order  to  determine  which  of  the  6  geometries  proposed 
is  the  best  design  to  minimize  the  entropy  produced. 

The  model  with  4  spirals  is  found  to  be  the  best  geometry  to  dis¬ 
tribute  the  gases.  The  main  advantages  of  this  model  are:  the  larger 
power  developed,  the  more  uniform  current  density  produced,  a 
uniform  management  of  water,  and  the  relatively  small  pressure 
drop.  The  model  with  8  spirals  produces  the  worst  performance  of 
the  fuel  cell.  The  rest  of  the  models  have  an  average  performance 
but,  still  lower  than  the  4-channel  performance. 

The  present  results  are  in  complete  disagreement  with  an  ear¬ 
lier  work  [14],  where  the  cell  performance  of  2  fuel  cells  with  1  and 
2  spiral’s  geometry  was  reported.  The  previous  result  established 
that  the  2-channel  (Fermat’s  spiral)  has  a  much  better  performance 
than  the  one  with  1  spiral.  The  present  results  show  a  similar  perfor¬ 
mance  between  the  two  models  and  there  is  not  an  improvement  of 
the  performance.  The  possible  differences  can  be  the  electrochemi¬ 
cal  parameters  used  or  the  boundary  conditions  that  are  completely 
different.  The  present  analysis  is  more  acceptable  due  to  the  fact 
that  the  inlet  gases  conditions  are  calculated  based  on  stoichio¬ 
metric  parameters. 

The  geometry  that  produces  the  least  amount  of  entropy  is  the 
1 -channel  fuel  cell  followed  by  the  2-channel  fuel  cell.  The  3  and  4- 
channel  fuel  cell  models  have  an  average  performance  as  compared 
to  the  other  four  geometries.  The  6  and  8-channels  models  are  the 
worst  geometries,  generating  the  maximum  amount  of  entropy. 

All  things  considered  for  the  present  problem  lead  to  the  follow¬ 
ing  conclusions: 


-  The  4-channel  geometry  is  the  best  option  to  implement  as  a  bipo¬ 
lar  plate.  It  generates  the  largest  current  density,  and  develops 
the  largest  power.  It  could  be  argued  that  it  does  not  produce  the 
lowest  amount  of  entropy  but  its  values  are  still  reasonable  good. 
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-  The  second  best  option  for  the  fuel  cell  geometry  is  the  3-channel. 
Its  performance  is  very  similar  to  that  of  the  4-channel  model  but 
with  the  difference  that  it  produces  a  lower  current  density. 

-  The  1 -channel  and  2-channel  models  have  an  average  perfor¬ 
mance  from  all  points  of  views,  being  good  options  if  the  criterion 
to  meet  is  the  minimum  entropy  generation. 

-  Finally,  the  6-channel  and  the  8-channel  develop  the  worst  per¬ 
formance  as  a  fuel  cell.  They  produce  the  lowest  current  density 
and  the  largest  entropy  production. 

Thus,  the  number  of  channels  used  considerably  affects  the  cell 
performance,  finding  the  optimal  value  between  3  or  4  channels.  It 
is  not  convenient  to  increase  the  number  of  channels  because  there 
is  not  a  significant  improvement  on  the  cell  performance,  and  also 
the  manufacturing  of  the  flow  channels  becomes  more  complex. 

Two  different  parameters  were  used  to  compare  the  cell  perfor¬ 
mance  based  on  the  entropy  analysis,  the  Bejan  number  and  the 
n  number.  The  present  results  are  in  complete  agreement  with 
those  previous  results  of  Damian  et  al.  [18]  in  that  the  matter  flow 
has  a  more  significant  effect  than  the  thermal  effects.  The  work  by 
Damian  and  the  presented  work  enhance  the  use  of  the  n  number 
as  the  best  option  to  analyze  and  characterize  fuel  cells. 
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